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1.  INTRODUCTION 


1.1  Purpose  of  Program 

This  report  describes  the  technical  objectives  and  reviews  the  research 
and  development  performed  under  Contract  DAAB07-74-C-0272.  The  project  began 
on  21  June  1974  and  has  been  sponsored  by  the  Radar  Technical  Area  of  the 
Combat  Surveillance  and  Target  Acquisition  Laboratory  within  the  U.  S.  Army 
Electronics  Command  (ECOM) . The  purpose  of  the  program  has  been  to  design 
and  build  an  Environment  and  Radar  Operation  Simulator  (EROS) . EROS  synthesizes 
controllable  radar  backscatter  combining  realistically  simulated  returns  from 
a variety  of  targets  and  clutter.  This  backscatter  incorporates  amplitude  and 
frequency  fluctuations  due  to  variations  in  radar  cross  section,  Doppler  shifts, 
and  other  phenomena.  The  resulting  synthetic  waveform  is  available  for  inser- 

tion into  the  radar  at  IF,  video,  and  Doppler. 

The  simulated  scenario  is  specified  by  the  EROS  operator  and  includes 
the  locations  and  movements  of  objects  producing  radar  backscatter:  troops, 
vehicles,  vegetation,  hills,  insects,  etc.  The  scenario  description  also 
includes  radar  and  reflectivity  characteristics:  antenna  pattern,  range 
dependence,  etc. 

The  EROS  system  translates  scenario  information  into  signal  data  in  a 
form  appropriate  for  real-time  simulation.  During  real-time  simulation,  the 
EROS  system  combines  scenario  data,  signals  from  the  subject  radar,  and  oper- 
ator commands  into  realistic  radar  backscatter,  which  is  applied  to  the  radar 
receiver.  This  form  of  simulation  permits  laboratory  testing  as  an  alternative 
to  field  testing  and  has  a number  of  potential  advantages: 

o EROS  testing  promises  to  be  less  expensive  than  field  testing, 
because  of  the  convenience  of  the  laboratory  as  a test  site. 

o The  results  are  controllable  to  a high  degree  of  scenario  detail. 
Unpredictable  influences  such  as  weather  changes  do  not  affect 
EROS  tests. 

® The  tests  are  repeatable.  This  has  obvious  benefits  in  experi- 
mentation and  in  repair  of  radar  receiver  components. 

° Signal  insertion  at  frequencies  other  than  RF  permits  simulated 
field  testing  of  specific  receiver  components  (e.g.,  embodying 
unproven  clutter  reduction  techniques)  without  having  to 
incorporate  this  component  in  a complete  radar. 

In  principle,  it  is  highly  desirable  for  EROS  to  be  applicable  to  all 
existing  troop  surveillance  radars  and  to  embody  sufficient  flexibility  to 
accommodate  as  yet  unspecified  radar  designs.  This  is  a long-range  goal  and 
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not  a specific  objective  of  the  current  program.  A particular  subject  radar, 
the  AN/PPS-15,  has  been  chosen  to  prove  feasibility  of  the  EROS  concept  and 
design,  and  EROS  has  been  built  to  address  this  radar.  However,  the  design 
incorporates  flexibility  that  will  permit  easy  adaptation  of  the  first  EROS 
system  to  other  radars. 


1.2  EROS  Description 

The  scenario  to  be  simulated  (see  Figure  1)  is  an  annular  shaped  region, 
which  represents  the  portion  of  the  battlefield  visible  to  the  radar.  The 
region  is  subdivided  into  concentric  rings  (at  constant  range  from  the  subject 
radar)  and  radial  columns  (at  constant  azimuth  with  reference  to  the  subject 
radar).  The  simulated  antenna  beam  consists  of  a contiguous  subset  of  m of 
these  azimuth  columns.  The  center  of  the  beam  is  determined  by  sampling  a 
voltage  from  the  subject  radar,  which  represents  its  antenna  angle. 

The  intersection  of  range  rings  and  azimuth  columns  defines  a sub- 
division of  the  scenario  into  four-sided  cells,  which  will  hereafter  be  called 
"simulation  cells."  The  concept  of  the  simulation  cell  is  useful  in  providing 
a reference  for  the  EROS  operator  in  defining  the  scenario.  Moreover,  much 
of  the  internal  communication  between  EROS  components  is  expressed  in  terms 
of  simulation  cells. 


The  EROS  hardware  consists  of  a digital  computer,  special  purpose 
digital  hardware,  and  analog  hardware.  The  flow  of  data  between  these  compo- 
nents is  illustrated  in  Figure  2.  The  computer  contains  a cathode  ray  tube 
graphic  display  unit  and  approximately  40  million  bits  of  online  disk  storage. 
These  facilities  are  used  to  assist  the  EROS  operator  in  defining  the  scenario 
and  to  participate  in  the  real-time  simulation  by  supplying  control  information 
to  the  digital  hardware.  The  digital  hardware,  in  turn,  generates  synthetic 
clutter  backscatter,  applies  antenna  pattern  weighting,  and  combines  clutter 
signals  with  target  signals.  Thn  combined  data  are  sent  to  the  analog  hard- 
ware. The  analog  hardware  applies  range-dependent  weighting  (e.g.,  R“^  power 
attenuation)  and  generates  a video  frequency  signal  incorporating  the  subject 
radar's  transmitter  modulations. 


The  target  and  clutter  signals  reflect  the  Doppler  and  radar  cross 
section  fluctuations  of  the  simulated  backscatter.  In  the  computer  and  digital 
hardware  these  are  expressed  in  terms  of  sampled  waveforms.  Each  sample  is  a 
complex  quantity  represented  by  the  two  quadrature  components  of  the  waveform. 

Target  signals  in  EROS  are  recorded  on  the  computer  disk  during  scenario 
definition.  Additional  target  data  (specified  via  the  computer's  graphic 
display  unit)  consists  of  target-movement  definition  and  a schedule  of  events. 
This  information  is  compiled  onto  a disk  file  and  is  subsequently  retrieved 
for  transmission  to  the  digital  hardware  during  real-time  simulation.  To 
ensure  realism,  the  target  signal  strength  is  apportioned  between  two  adjacent 
range  rings.  More  precisely,  assume  that  the  target's  signal  strength  is  s 
and  that  the  target  is  between  the  center  of  range  ring  r and  the  center  of 
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Figure  2. 


EROS  Data  and  Control  Flow. 
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range  ring  r + 1.  Let  dQ  and  d^^  denote  the  respective  distances  between 
target  and  the  centers  of  range  rings  r and  r + 1.  Then  the  backscatter 
the  target  is  represented  as 

dj^s  ^0^ 

— r ; — r~  from  range  ring  r and  — ; ; — r-  from  range  ring  r + 1. 

do  + d^  do  + dj^ 

The  reason  for  doing  this  is  to  avoid  an  abrupt  discontinuity  (in  range-gated 
radar  reception)  as  the  target  crosses  the  boundary  between  range  rings.  A 
similar  problem  in  azimuth  boundaries  is  avoided  by  utilizing  very  high 
resolution  in  the  representation  of  the  antenna  pattern. 


Clutter  is  synthesized  in  the  digital  hardware  by  applying  pseudo-random 
sequences  (digital  noise)  to  digital  filters.  The  outputs  from  the  filters  have 
spectra  and  amplitude  distributions  that  are  controllable  as  functions  of  the 
filter  parameters.  Clutter  backscatter  is  represented  as  the  sum  of  the  radar 
returns  from  moving  scatterers  (e.g.,  foliage,  insects,  etc.)  and  fixed  scat- 
terers  (e.g.,  ground,  buildings,  boulders,  etc.).  The  return  from  the  moving 
portion  of  clutter  is  created  by  two  identical  2-pole  filters,  which  generate 
the  two  quadrature  components.  The  return  from  the  fixed  portion  of  clutter 
is  simply  a complex  constant.  The  2-pole  filter  employs  two  feedback  multi- 
pliers and  02.  The  state  of  the  filter  is  defined  by  the  most  recent  output 
w^  and  the  previous  output  W2.  The  transition  to  the  next  filter  state  is 
described  by  the  computation: 


w'  ■<'  w;j^  + 02  W2  + X 
W2  w^ 

W]^  w' 

where  w'  is  the  next  filter  output,  and  x is  the  pseudo-random  input. 


(1-1) 


The  computations  are  performed  by  two  identical  units  of  hardware  for 
the  two  quadrature  components.  The  parameters  and  state  information  are  main- 
tained in  read-write  memory — the  memory  holds  one  set  of  such  data  for  each 
simulation  cell  in  the  antenna  beam.  The  filter  hardware  performs  computations 
in  an  endless  loop: 


update  the  state  of  filter  number  1 

update  the  state  of  filter  number  2 

update  the  state  of  filter  number  k (k  denotes 

the  number  of  filters  in  the  beam) 

update  the  state  of  filter  number  1 
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This  arrangement  avoids  maintaining  k sets  of  hardware  for  the  k simulation 
cells.  Moreover,  by  permitting  the  digital  computer  to  write  into  the  para- 
meter memory,  it  provides  a means  for  simulating  the  clutter  response  to 
antenna  scanning.  As  simulation  cells  drop  off  the  trailing  edge  of  the 
antenna  beam,  new  cells  enter  the  leading  edge.  The  computer  overwrites  the 
parameters  of  the  trailing-edge  cells  with  the  parameters  of  the  leading-edge 
cells. 


One  loop  through  the  k cells  constitutes  one  Doppler  sampling  cycle. 
The  clutter  hardwar"  produces  samples  for  k cells  during  the  Doppler  sampling 
cycle.  These  are  combined  with  a small  number  (at  most  two)  of  target-signal 
samples  with  their  respective  antenna-pattern  weights  and  range  ring  numbers 
that  the  computer  has  sent.  These  data  are  combined  as  follows: 

° Azimuth  weights  are  applied  to  the  target  and  clutter  samples 

(the  antenna  pattern  for  clutter  is  sent  from  the  computer  during 
initialization  and  the  digital  hardware  stores  this  array  in  a 
random-access  memory). 

o All  clutter  samples  at  the  same  range  are  added  together  (azimuth 
integration) . 

° The  target  samples  are  added  to  the  appropriate  clutter  sums 
(based  on  the  specified  range-ring  number). 


The  result  is  range  profile — one  value  for  each  range  ring. 

The  range  profile  is  recomputed  at  the  Doppler  sampling  frequency,  which 
is  several  times  slower  than  the  pulse-repetition  frequency  (PRF)  of  the 
AN/PPS-15  radar.  Therefore,  the  digital  hardware  stores  each  range  profile 
in  a random-access  memory.  This  same  data  is  retrieved  at  the  PRF;  several 
repetitions  of  the  same  data  are  used  before  a new  range  profile  is  stored. 


The  range  profile  is  sent  to  the  analog  hardware  at  the  PRF.  Each  trans- 
mission is  delayed  a fixed  interval  after  the  radar's  pulse  trigger;  this  delay 
corresponds  to  the  distance  between  the  radar  and  the  closest  simulated  range. 
The  range  samples  are  sent  in  ascending  range  order.  The  Interval  between 
.•mples  corresponds  to  the  width  of  the  simulated  range  ring.  Each  range 
sample  is  accompanied  by  the  corresponding  range-dependent  attenuation.  For 
example,  the  range  attenuation  might  vary  proportionally  to  R~^  to  simulate 

power  loss.  After  the  transmission  of  the  range  sample  for  the  most  distant 
ring,  the  digital  hardware  waits  for  the  trigger  denoting  the  next  radar  pulse. 


The  analog  hardware  employs  3 D/A  converters:  1 for  the  in-phase  signal 
component,  1 for  the  quadrature  signal  component,  and  1 for  the  range  weight. 

The  range  weight  is  applied  to  the  two  signal  components  using  an  analog  multi- 
plier. The  signals  are  then  fed  through  an  analog  integrator  in  order  to 
simulate  the  transmitter  modulation  of  the  AN/PPS-15  radar.  Finally,  a 
controllable  amount  of  analog  noise  is  added  to  account  for  antenna  and  receiver- 
front-end  noise,  and  the  result  is  applied  to  the  video  stage  of  the  AN/PPS-15 
radar. 
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2. 


EROS  PARAMETERS 


This  section  summarizes  decisions  that  were  made  concerning  parameters 
for  the  first  EROS  implementation  model.  Many  of  these  parameters  are  based 
on  characteristics  of  the  AN/PPS-15  radar,  and  they  should  be  reviewed  care- 
fully when  considering  EROS  revisions  to  address  a different  radar. 


2.1  Video  Sampling  Rate — Range  Resolution 

The  video  sampling  rate  refers  to  the  transmission  of  range-ring  samples 
from  the  EROS  digital  hardware  to  the  EROS  analog  hardware.  This  rate  must  be 
fast  enough  to  yield  a range  of  resolution  of  approximately  15  meters,  one-half 
the  range  resolution  of  the  AN/PPS-15  radar.*  The  sampling  rate  is  also  chosen 
for  convenience  to  be  a multiple  of  the  37.5  KHz  rate  of  the  AN/PPS-15  trans- 
mitter's modulator  trigger.  The  value  selected  equals  (266)  x (37.5  KHz)  or 
9.975  MHz.  In  summary: 


Video  sampling  rate 
Video  sampling  period 
Range  ring  width** 


= 9.975  MHz 

= 100.25  ns 

= 15.027  meters. 


2.2  Doppler  Sampling  Rate 

The  AN/PPS-15  radar  has  a Doppler  cut-off  frequency  of  about  1 KHz. 
Therefore,  a sampling  rate  on  the  order  of  4 KHz  is  judged  to  be  adequate  for 
faithful  representation  of  target  and  clutter  envelopes.  Since  it  is  conve- 
nient to  choose  the  sampling  rate  to  be  an  Integral  divisor  of  the  37.5  KHz 
modulator  trigger  rate,  the  sampling  frequency  (■^)  x (37.5  KHz)  = 4.167  KHz 
has  been  selected.  In  summary: 

Doppler  sampling  rate  = (^)  x (37.5  KHz)  = 4.167  KHz 

T = Doppler  sampling  period  = 240  pseconds. 


*A  sampling  rate  corresponding  to  the  30  meter  range  resolution  is  Inadequate 
because  of  the  representation  of  targets  in  two  adjacent  rings. 

**This  value  is  based  on  the  radar  propagation  speed  of  2.99796  x 10^ 
meters/sec . 
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2.3  Number  of  Range  Rings 

The  AN/PPS-15  is  sensitive  to  radar  return  from  ranges  between  50  meters 
and  3000  meters.  Simulating  this  2950  meter  range  would  require  190  15-meter 
rings.  However,  it  is  not  necessary  to  simulate  the  entire  range  in  any  one 
simulation  run.  EROS  will  provide  the  capability  of  simulating  64  contiguous 
rings  during  any  one  simulation  run,  and  the  specific  64  rings  are  selectable 
as  a simulation  parameter.  In  summary: 

The  number  of  range  rings  = 64 

The  total  simulated  range  = (64)  x (15.027)  = 962  meters. 

2.4  Azimuth  Width  of  Simulated  Clutter 


The  azimuth  antenna  pattern  for  the  AN/PPS-15  indicates  that  the  null 
following  the  first  sidelobe  occurs  approximately  11°  from  center  beam.  There- 
fore, a simulated  area  22°  wide  would  include  the  main  beam  and  the  first 
sidelobe  on  both  edges.  A convenient  width  for  the  EROS  simulation  is  of 
a circle  = 22.5°  = 400  mils.  In  summary: 

Azimuth  width  of  simulated  clutter  = 22.5°  = 400  mils. 


2.5  Azimuth  Width  of  Clutter  Cells 

The  azimuth  column  width  in  the  scenario  partition  determines  the 
azimuth  width  of  each  simulated  clutter  cell.  A convenient  value  for  the 
azimuth  column  width  is  of  a circle  = 25  mils  = 1.4°.  Since  the  azimuth 
antenna  pattern  for  the  AN/PPS-15  indicates  that  the  main  beam  is  about  6°  wide 
between  the  3 dB  points  and  about  13°  wide  between  the  nulls,  this  would  provide 
at  least  4 azimuth  columns  between  the  3 dB  points  and  at  least  9 azimuth  col- 
umns between  the  main-beam  nulls.  This  is  judged  to  be  adequate  for  realistic 
clutter  simulation.  In  summary: 

Azimuth  width  of  clutter  cell  = 25  mils  = 1.4°. 


2.6  Range  Width  of  Clutter  Cells 

The  range  width  of  clutter  cells  has  been  chosen  to  be  2 range  rings 
or  approximately  30  meters.  This  provides  32  clutter  cells  for  each  azimuth 
column.  Subdividing  the  simulated  clutter  scenario  into  32  30-meter  rings  is 
judged  to  provide  adequate  realism  and  flexibility.  In  summary: 

Range  width  of  clutter  cell  = (2)  x (15.027)  = 30.054  meters 

Number  of  clutter  cells/azimuth  column  = 32. 
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2.7  Azimuth  Dependence  Upon  Target  Samples 

The  antenna  pattern  weights  to  be  applied  to  target  samples  are  stored 
in  a table  in  the  computer's  main  memory.  Because  of  real-time  constraints, 
this  table  is  not  interpolated;  instead,  the  entry  that  corresponds  to  the 
closest  target  azimuth  is  used.  If  the  azimuth  spacing  is  too  coarse,  the 
antenna  gain  will  remain  constant  for  several  Doppler  sampling  cycles  and  then 
undergo  a jump  discontinuity.  Repetition  of  this  behavior  at  low  frequencies 
will  introduce  an  annoying  spurious  tone  in  the  aural  output  of  the  AN/PPS-15 
radar.  A study  was  performed  to  determine  how  frequently  the  antenna  pattern 
gain  must  change  in  order  to  avoid  this  spurious  signal.  The  result  is  that 
the  antenna  gain  must  change  at  least  1244  times/second  to  reduce  the  worst- 
case  spurious  signal  to  72  dB  below  the  saturating  signal  (details  are  supplied 
in  the  Appendix).  This  constraint  coupled  with  the  90  mils/second  scan  rate 
of  the  AN/PPS-15  radar  implies  an  azimuth  spacing  of  2“^^  circles.  The  width 
of  the  antenna  pattern  is  of  a circle,  and  pattern  symmetry  implies  that 
only  half  of  the  pattern  needs  to  be  repeated  in  the  table.  In  summary: 

Azimuth  spacing  of  antenna-pattern  table  = 2 

circles  = .00275° 

Number  of  points  in  pattern  table  (-^  pattern)  = 2^^  = 4096 


2.8  Sample  Representation 

Target  and  clutter  signals  are  represented  by  bit  patterns  in  the 
computer  and  in  the  digital  hardware.  Two's  complement  representation  is 
employed,  and  the  quantities  are  generally  interpreted  as  binary  fractions. 

Thus  the  bit  pattern  bQbj^  . . . bj^  (bj^  equals  0 or  1)  is  interpreted  as  the 
value 

-bo  + (^)  b^  + C-^)  b^  + . . . + (-j)™  b^. 

The  interface  between  the  digital  hardware  and  the  analog  hardware 
employs  12-blt  D/A  converters.  This  permits  72  dB  dynamic  range  in  the  signal 
amplitude  due  to  Radar  Cross  Section  (RCS)  variation.  (Note  that  the  addition- 
al dynamic  range  due  to  range  dependence  is  not  constrained  by  this  choice, 
because  range  weights  are  transmitted  through  a separate  D/A  converter.) 

The  RCS  0 of  a signal  sample  is  proportional  to  the  squared  amplitude 
of  the  sample.  The  values  I and  Q of  the  in-phase  and  quadrature  components 
of  the  sample  are  determined  by  using  the  binary-fraction  interpretation 
stated  above,  and  the  proportionality  constant  is  denoted  by  K^.  Thus, 

o = (I^  + q2)  (2-1) 

2 

The  constant  20,000  m has  been  chosen  for  K^.  In  summary. 
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Bit  patterns  are  interpreted  as  2's  complement  binary 
fractions. 

Twelve-bit  D/A  converters  are  used  for  each  quadrature 
component. 

The  proportionality  between  squared  amplitude  and  RCS  is 

20,000  m2. 


2.9  Antenna  Azimuth  Sampling  Rate 

The  antenna  azimuth  accuracy  of  the  AN/PPS-15  radar  is  judged  to  be 
roughly  the  same  as  the  read-out  resolution,  5 mils.  The  scan  rate  of  the 
AN/PPS-15  radar  is  90  mils/sec.  EROS  selects  samples  from  the  antenna 
azimuth  voltage  at  about  the  same  rate  as  the  frequency  at  which  the 
AN/PPS-15  azimuth  changes  by  the  accuracy  amount,  90/5  = 18  times/sec. 

Between  samples  the  azimuth  value  is  computed  by  linear  interpolation.  If 
the  azimuth  were  sampled  less  frequently,  the  azimuth  accuracy  of  EROS  would 
be  degraded.  If  the  azimuth  were  sampled  more  frequently,  the  discontinuities 
due  to  accuracy  limitations  in  the  azimuth  voltage  would  be  more  pronounced. 
The  possible  adverse  consequences  of  the  latter  problem  are  spurious  signals 
at  the  azimuth  sampling  rate,  an  effect  similar  to  that  described  in 
Section  2.7  due  to  discontinuities  in  the  antenna-pattern  table.  The  azimuth 
sampling  frequency  has  been  chosen  for  convenience  to  be  1/256  the  Doppler 
sampling  frequency  or  about  16.28  Hertz.  In  summary 

, ^ ....  1 . ^ Doppler  Sampling  Rate 

Antenna  azimuth  sampling  rate  = — “ 

/3D 

= 16.28  Hertz 
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3.  CLUTTER  SIMULATION 


Backscatter  Representation 


Radar  backscatter  is  a modulated  RF  signal.  Radar  cross  section  fluc- 
tuations introduce  amplitude  modulations,  and  the  Doppler  shifts  due  to  move- 
ments of  reflecting  elements  constitute  frequency  modulations.  In  general, 
the  bandwidth  of  this  modulation  is  very  narrow  compared  to  the  transmitted 
RF.  This  narrow  bandwidth  permits  a much  simpler  representation  of  the 
backscatter  in  terms  of  its  complex  envelope.* 

The  complex  envelope,  expressed  in  terms  of  the  in-phase  component  and 
quadrature  component,  is  the  radar  backscatter  representation  used  in  EROS. 
Since  the  video  signal  of  the  AN/PPS-15  radar  is  produced  by  a coherent 
detector,  the  complex  envelope  can  be  applied  directly  to  the  video  stage 
of  the  radar.  For  radars  employing  non-coherent  detectors  it  would  not  be 
correct  to  apply  the  complex  envelope  at  video  without  additional  EROS  pro- 
cessing. Instead  it  might  be  necessary  to  create  a modulated  signal  from 
the  synthetic  envelope  and  either  (1)  to  insert  this  as  RF  or  IF,  or  (2)  to 
replicate  the  radar's  detector  in  order  to  insert  the  detected  signal  into 
the  radar. 

To  be  realistic,  the  envelope  of  the  radar  backscatter  should  be  re- 
presented as  a random  process.  However,  since  EROS  is  designed  to  produce 
repeatable  results,  the  simulated  backscatter  envelope  cannot  be  truly 
random.  Nevertheless,  EROS  does  produce  a facsimile  of  a typical  waveform 
selected  from  the  random  process  specified  by  the  scenario  parameters. 

The  characteristics  of  a random  process  are  often  expressed  statis- 
tically. Two  commonly  cited  characteristics  are  the  amplitude  distribution 
and  the  power  spectral  density.  For  a stationary  continuous-valued  random 
process  x(t),  the  amplitude  distribution  is  the  probability  distribution  of 
|x(tQ)|,  where  tQ  is  any  fixed  value  of  time  t.  The  power  spectral  density 
of  x(t)  is  the  Fourier  transform  of  the  auto-correlation  function  of  x(t). 
(Refer  to  [1]).  The  clutter  model  in  EROS  is  a random  process  with  controllable 
amplitude  distribution  and  power  spectral  density.  The  underlying  assumption 
of  this  model — that  backscatter  can  be  represented  as  a stationary  random 
process — is  not  necessarily  correct.  There  is  a risk  that  the  use  of  sta- 
tionary statistics  to  describe  clutter  will  detract  from  the  realism  of  the 
simulation.  This  question,  among  others,  is  necessarily  deferred  to  an 
evaluation  of  EROS. 


*Section  4 of  the  First  EROS  Quarterly  Report  gives  a more  detailed 
discussion. 
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EROS  hardware  produces  synthetic  backscatter  in  digital  form.  This 
is  satisfactory  as  long  as  the  sampling  frequency  is  several  times  the  band- 
width of  the  random  process.  The  carrier  frequency  of  the  AN/PPS-15  radar  is 
about  10  GHz.  Since  the  bandwidth  of  X-band  clutter  is  at  most  hundreds  of 
Hertz,  and  since  EROS  simulation  will  restrict  target  Doppler  frequencies 
to  less  than  1 KHz,  the  4.167  KHz  Doppler  sampling  frequency  is  adequate. 
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3.2  Clutter  Amplitude  Distribution 

Of  the  various  models  available  that  describe  clutter  amplitude  dis- 
tributions, the  Ricean  has  been  chosen  for  EROS  implementation.  This  decision 
was  prompted  by  several  factors: 

° The  Ricean  density  function  is  based  on  a physical  model  of 
ground  clutter. 

° The  Ricean  model  nas  widespread  support  in  the  literature  as 

a realistic  amplitude  distribution  model  for  ground  clutter.  It 
has  been  successfully  curve-fitted  to  clutter  measurement  data 
taken  under  a variety  of  environmental  conditions. 

° The  assumption  of  a Ricean  amplitude  distribution  makes  the 

clutter  synthesis  implementation  manageable.  The  SAM-D  project 
incorporates  a Ricean  clutter  model  in  its  radar  simulator  [2]. 


A derivation  of  the  Ricean  distribution  from  a physical  model  of 
clutter  is  provided  by  Goldstein  [3]  in  which  a signal  is  formed  by  adding 
two  signal  components.  The  complex  envelope  of  one  of  the  signal  components 
has  a constant  magnitude  |s|.  The  complex  envelope  of  the  other  signal  com- 
ponent fluctuates  with  a uniform  phase  distribution  and  a Rayleigh  distribu- 
tion of  the  amplitude  |r|: 

2|r|  / |r|2  \ 

) = exp  I I (3- 

Pq  \ Pq  / 

where  denotes  a probability  density  function. 


The  technique  for  generating  a random  process  with  a Ricean  amplitude 
distribution  is  based  on  Goldstein's  physical  model:  the  sum  of  two  signal 
components,  one  of  which  has  a constant  complex  envelope  S,  and  the  other, 
which  is  a random  process  with  complex  envelope  R(t),  is  described  as  follows. 
The  random  process  R(l)  is  stationary  with  complex  values  of  the  form 
X(t)  + jY(t).  The  probability  distribution  of  RCCq)  for  any  fixed  time  tQ 
is  a bivariate  Gaussian  distribution  of  the  random  phasor  X(to)  + j^Cto;  with 
density  function 


Wi(X,Y) 


exp 


X^  + y2  \ . 

PO  I 


(3-2) 
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The  implementation  parallels  this  physical  model.  Two  series  of  inde- 
pendent random  samples  . . . uq,  U2>  • • • and  . . . v_]^,  vq,  Vj^,  V2, 

. . . are  fed  through  two  identical  real  discrete  filters  resulting  in  two 
real  discrete  random  processes  Yj(nT)  and  Yg(nT),  respectively.  These  are  then 
added  to  the  complex  constant  A to  form  the  discrete  complex  envelope  H(nT)  of 
clutter : 

N(nT)  = A + Yj.  (nT)  + jYg  (nT)  . (3-3) 

It  is  difficult,  in  general,  to  determine  the  probability  distribution 
of  Yj  and  Yq.  In  the  particular  case  when  u^  and  v^  are  Gaussian  samples, 
it  can  be  proven  that  Yj  and  Yq  are  also  Gaussian  distributed.*  Results  show 
that  these  filter  output  statistics  remain  Gaussian  under  considerably  wider 
circumstances.  Experiments  with  filters  typical  of  those  used  in  EROS  for 
clutter  synthesis  have  indicated  that  the  probability  distribution  is  approxi- 
mately Gaussian  even  if  the  inputs  are  randomly  selected  from  the  set  (-1,  1}. 

A pseudo-random  number  generator  with  a two-point  distribution  is  much  easier 
to  implement  than  one  whose  distribution  is  pseudo-Gaussian.  Therefore,  the 
two-point  random  number  generator  has  been  adopted  for  EROS  implementation. 

The  mean  of  the  distribution  of  Yj  and  Yq  is  not  zero  in  general 
(clearly  Yj.  and  Yy  are  identically  distributed  and  therefore  have  the  same 
mean).  This  is  true,  in  spite  of  the  fact  that  the  pseudo  random  inputs  have 
zero  mean,  because  of  truncation  errors  introduced  by  the  filter's  multipliers. 
However,  this  does  not  increase  the  complexity  of  the  filter  hardware,  since 
the  non-zero  bias  can  be  compensated  by  the  appropriate  choice  of  A.  Let  u 
denote  the  mean  of  Yj.  and  Yq.  Then  (3-3)  can  be  rewritten  in  the  form 

N(nT)  = A + m(1  + j)  + [Yj(nT)  - tj]  + j [YQ(nT)  - p].  (3-4) 

Yj-(nT)  - p and  YQ(nT)  - p are  independent  and  identically  distributed  Gaussian 

random  variable  with  zero  mean.  Their  bivariate  density  function  has  the  form 

(3-2).  Therefore,  iN(nT) j has  a Ricean  distribution. 

The  term  A + p(l  + j)  represents  the  component  of  clutter  backscatter 
due  to  fixed  scatterers.  Its  amplitude  should  correspond  to  the  radar  cross 
section  0^  of  this  component  of  clutter.  From  (2-1) 

°1  = '^o  + ‘'(1  + j)  . (3-5) 


*This  is  because  Yj  and  Yq  can  be  represented  as  convolutions  of  the  discrete 
filter's  unit-impulse  response  with  a random  sequence.  This  convolution  is, 
in  turn,  simply  a linear  combination  of  the  random  inputs.  A linear 
combination  of  Gaussian  random  variables  yields  a Gaussian  random  variable. 


The  term  [Yj  (nT)  - uj  +j[YQ  (nT)  - M]  represents  the  component  of  clutter 
backscatter  due  to  moving  scatterers.  Its  amplitude  should  correspond  to  the 
radar  cross  section  '^2  moving  part  of  clutter.  However,  since  the 

radar  cross  section  from  a collection  of  moving  scatterers  is  not  constant, 

02  more  appropriately  denotes  the  average  radar  cross  section. 


= Kg  average 


{[Yj  (nT)  - 


+ 


[Yq  (nT) 


(3-6) 


= Kg  variance  (Y^.  (nT)  + jY^  (nT) } 


3.3  Clutter  Spectral  Distribution 

The  clutter  spectral  model  adopted  for  EROS  is  expressed  by  the 
formula* 

P(f)  = — , (3-7) 


where  P(f)  denotes  the  power  spectral  density  as  a function  of  frequency, 
fg  is  the  cutoff  frequency,  and  n expresses  the  rate  of  high  frequency  roll- 
off. The  value  of  n has  been  observed  to  vary  between  2 and  4 depending  upon 
the  transmitter  frequency  and  the  type  of  clutter  producing  backscatter  [ 5] . 
Formula  (3-7)  applies  only  to  the  fluctuating  component  of  clutter.  This 
model  was  derived  by  fitting  curves  to  measurement  data  from  which  the  0-fre- 
quency component  was  filtered  out.  Note  also  that  (3-7)  is  normalized  so  that 
radar-cross  section  dependence  is  eliminated,  and  P(f)  approaches  1 as  f 
approaches  0.  Thus  the  spectral  requirement  for  the  EROS  clutter  model  is 
that  the  power  spectral  density  of 

[Yi(nT)  - p]  + j[Y^;j(nT)  - p] 

be  proportional  to  P(f)  as  expressed  in  (3-7).  The  following  sections  of  this 
report  describe  how  this  has  been  implemented. 


3.4  Filter  Design 

It  is  well  known  that  analog  filters  can  be  used  to  shape  white  noise 
into  a signal  with  a desired  spectrum.  A similar  property  applies  to  digital 
filters.  Let  Pjj(u))  denote  the  power  spectral  density  of  a random  process  x 
as  a function  of  frequency  u (in  units  of  radians/second) 

0)  = 2Trf  . (3-8) 


*This  model  is  suggested  by  [4]. 


! 


1 


\ 

*i 


i 

? 

f 

i 

9 

I 

1 


j 

5 


I 


If  a random  process  x is  applied  to  a digital  filter  whose  unit-impulse  response 
has  the  z-transform  H(z),  then  the  output  random  process  y has  the  power  spectral 
density  Py(m)  given  by* 

P^(w)  = |H(e^”‘')  1^  P (w)  . (3-9) 

Equation  (3-9)  expresses  the  motivation  for  using  digital  filters  to  synthesize 
the  desired  spectrum.  If  the  input  random  process  x has  an  autocorrelation 
function  with  a constant  value  at  non-O-lag  and  with  a different  value  at  0-lag, 
then  the  spectrum  of  the  input  process  is  constant  except  for  integer  multiples 
of  the  sampling  frequency.  Therefore,  the  spectrum  of  the  output  process  y is 
proportional  to  |H(ej''^)|2  except  at  multiples  of  the  sampling  frequency. 

The  sampling  frequency  is  sufficiently  high,  so  that  Che  only  multiple  within 
the  pass  band  of  the  subject  radar  is  the  0-frequency  term.  This  term  repre- 
sents part  of  the  backscatter  due  to  stationary  scatterers:  the  constant 
term  to  be  added  to  the  filter  output  accounts  for  tne  rest  of  this  backscatter. 

The  approach  adopted  for  the  clutter  filter  design  has  been  to  select 
an  analog  filter  that  would  provide  the  specified  spectral  shape  (3-7)  and  to 
derive  the  digital  filter  from  the  analog  filter.  Consider  the  transfer  func- 
tion H^(s)  of  an  analog  filter.  The  spectral  response  of  an  analog  filter  is 
obtained  by  evaluating  the  squared  magnitude  of  H3(s)  along  the  imaginary  axis. 
Therefore,  the  requirement  that  the  analog  filter  have  the  specified  spectral 
shape  is  restated  in  the  form 


|H  (jm)!^  = — (3-10) 

a n 

c 

where  to  is  given  by  (3-8)  and  the  proportionality  constant  b denotes  the 
singal  strength  of  the  fluctuating  component  of  clutter  backscatter.  The  cut- 
off frequency  in  radians  per  second  is  denoted  by 

0)  = iTTf  (3-11) 

c c 


*Refer  to  Section  5.3  of  the  Second  EROS  Quarterly  Report.  To  be  perfectly 
accurate,  the  derivation  in  this  report  is  inexact  even  though  the  conclusion. 
Formula  (3-9),  is  still  correct.  Implicit  in  this  derivation  was  the  assump- 
tion that  the  z-transform  of  the  autocorrelation  function  converges  at  all 
points  on  the  unit  circle.  This  assumption  is  not  valid,  because  the  input 
sequence  has  a periodic  component  in  the  form  of  a non-zero  mean  due  to 
multiplier  truncation  errors  (see  Section  2.3  of  the  Sixth  EROS  Quarterly 
Report).  However,  if  the  autocorrelation  function  is  separated  into  periodic 
and  aperiodic  componenets,  and  if  the  Fourier  transform  is  applied  directly 
to  the  periodic  component  without  employing  the  z-transform,  the  result  (3-9) 
is  obtained. 
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Equation  (3-10)  suggests  a transfer  function  of  the  form: 


H (s)  = 


b s^S2 

(s  - Sj^Xs  - S2) 


(3-13) 


The  reasons  for  this  choice  are  as  follows: 


If  S]^  = cocexp(-jl35°)  — a 2 
satisfied  exactly  for  the  case  n = 
aliasing  is  expected  in  the  digital 
would  be  to  reduce  the  rate  of  high 
an  approximation  for  the  case  n=  3. 
mentation  a two-pole  filter  can  be 
one  of  the  feedback  multipliers  to 
solution  for  the  case  n = 2, 


-pole  Butterworth  filter  — then  (3-10)  is 
4.  Moreover,  a small  amount  of  frequency 
filter  derived  from  (3-13).  The  effect 
frequency  rolloff,  and  it  might  provide 
Finally  for  some  digital  filter  imple- 
converted  to  a one-pole  filter  by  setting 
zero.  A one-pole  filter  provides  an  exact 


In  order  for  (3-13)  to  correspond  to  a real-valued  unit-impulse  res- 
ponse, sj^  and  S2  must  either  be  complex  conjugates  or  real.  Moreover,  the 
filter  must  be  stable  which  implies  that  and  S2  be  in  the  left  half  of 
the  complex  plane. 

Applying  the  impulse- invarient  method  [6]  to  the  transfer  function 
(3-12)  yields  a digital  filter  whose  z-transform  has  the  form 


H(z)  = 


(3-14) 


ct^  z - a2Z 


The  parameter  3 corresponds  to  the  signal  strength  and  is  determined  by  the 
radar  cross  section  of  the  fluctuating  component.  The  parameters  and  02 
influence  the  spectral  behavior  of  the  clutter,  since  its  simulated  power 
spectral  density  (by  (3-9))  will  be  proportional  to  the  squared  magnitude 
of  H(z)  evaluated  on  the  unit  circle.  To  perform  such  evaluations,  it  is 
useful  to  express  and  02  in  terms  of  the  cutoff  frequency  the  pole 
angle  0,  and  the  sampling  Interval  T. 

If  the  analog  poles  are  complex  conjugates 

s.  = u)  e-^ 

1 P 


s„  = u)  e 
2 P 
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then  the  impulse  invarient  derivation  yields 


„ (Ur>T  COS6  / -r  . Q\ 

= 2e  P cos(WpT  sin  9) 


(3-15) 


2u»d  T cos  6 
= -e  ^ 


(3-16) 


The  magnitude  of  the  pole  is,  in  turn,  related  to  by  applying  (3-12): 


0)  = (1)  /cos  26  + /].  + cos^  26 

P c 


(3-17) 


On  the  other  hand,  if  the  poles  are  real 


s,  = - 0)  + D 

1 P 


s = - 0)  - D. 

2 P 


then  the  impulse  Invarient  derivation  yields 


= 2e““P^  cosh  DT 
-2m  T 

02  = -e  P • 


Again  by  applying  (3-11) : 


u)  = + 0)  ^ + 1/2  w ^ + 4 w 

p c c c 


(3-18) 

(3-iy) 


(3-20) 


A graphic  comparison  of  the  filter's  power  spectral  density  with  (3-7) 
was  performed  with  the  support  of  the  EROS  computer's  graphic  display  unit  and 
light  pen.  It  was  observed  that  pole  angles  between  152°  and  160°  provide  a 
good  approximation  to  (3-7)  for  the  case  n = 3.  The  high  frequency  rolloff  is, 
as  expected,  closer  to  the  desired  behavior  due  to  the  frequency  aliasing. 

In  order  for  a digital  filter  to  be  stable,  the  poles  must  be  within  the 
unit  circle.  Stability  is  important  in  this  application  to  prevent  the  feed- 
back values  from  overflowing.  The  stability  condition  applied  to  (3-14)  yields 
the  following  inequalities: 


02  > -1 

01  + 02  < 1 

02  - Oi  < 1 


(3-21) 


I 
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A more  detailed  analysis  of  filter  behavior  for  various  choices  of 
and  02  given  in  Section  2.2  of  the  Sixth  EROS  Quarterly  Report.  Included 
is  a map  of  the  (aj^,a2)  plane  indicating  correspondence  between  certain  points 
in  this  plane  with  parameters  of  the  analog  filter  (3-13)  from  which  (3-14) 
was  derived.  Other  points  in  the  (02^,02)  plane  are  identified  with  1-pole 
filter  behavior,  which  simulates  (3-7)  for  the  case  n = 2. 

The  network  corresponding  to  (3-14)  is  shown  in  Figure  3.  The  stabil- 
ity conditions  (3-21)  imply  that  [021  <1  and  that  lo-j^l  < 2.  Furthermore, 
(3-15)  suggests  that  negative  values  of  o^  will  not  be  particularly  useful. 
Therefore,  the  additional  restriction 

0 < ci^  < 2 (3-22) 

has  been  imposed.  This  permits  a convenient  modification  of  the  filter  to  the 
network  shown  in  Figure  4,  which  is  easier  to  implement  than  the  network  of 
Figure  3.  Both  multipliers  in  Figure  4 can  use  the  same  number  representation 
(the  representation  specified  in  Section  2.8),  since  the  factors  are  restrict- 
ed to  magnitudes  less  than  1.  Sixteen-bit  precision  has  been  selected  to 
obtain  sufficient  resolution  in  representation  of  of  a2  and  to  avoid  over- 
flow. 


3.5  Clutter  Azimuth  Weighting 

The  angular  width  of  the  antenna  beam  is  identical  to  the  width  of  16 
adjacent  clutter  columns  (see  Sections  2.4  and  2.5).  The  position  of  each 
clutter  column  relative  to  the  antenna  beam  varies  as  the  antenna  scans.  An 
illustration  of  beam  scanning  to  the  left  is  contained  in  Figure  5.  The  heavy 
curve  indicates  the  antenna  pattern  relative  to  clutter  cell  positions  when 
the  pattern  endpoints  are  exactly  aligned  to  span  16  columns.  The  seventeenth 
column  at  the  left  end  contains  new  filter  parameters  that  have  very  recently 
been  supplied  from  the  EROS  computer.  At  the  time  represented  by  the  heavy 
curve,  the  32  filters  in  this  seventeenth  column  are  generating  samples; 
however,  the  antenna  pattern  above  the  seventeenth  column  is  zero,  which 
nullifies  their  outputs.  The  dashed  curve  in  Figure  5 denotes  the  antenna 
pattern  shortly  thereafter.  At  this  time  the  antenna  pattern  weighting  has 
changed  for  all  of  the  cells,  and  the  weighting  for  cells  in  the  seventeenth 
column  is  now  nonzero. 

The  number  of  clutter  cells  simulated  at  any  one  time  equals 
(17)x(32)  = 544.  Each  of  the  544  cells  will  be  sampled  and  azimuth  weighted 
at  the  4.167  KHz  Doppler  sampling  rate,  i.e.,  within  a period  240  ysec.  These 
544  weighted  samples  are  subdivided  into  32  sets  of  17  samples;  each  set  of 
17  samples  is  taken  from  17  cells  at  the  same  range.  The  17  weighted  samples 
are  added  together  (azimuth  integration) , and  the  sum  is  the  total  clutter 
component  of  the  simulated  backseat  ter  from  that  range. 

The  clutter  azimuth  weights  are  maintained  in  a table  stored  in  digital- 
hardware  random-access  memory.  The  angular  spacing  between  points  in  this 
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table  is  2 circles  (roughly  the  same  as  the  5-mil  resolution  of  the 
AN/PPS-15's  antenna  read-out).  The  azimuth  weight  at  each  of  these  points 
is  given  by  the  formula 


w(e^)  = 


(3-23) 


where 


g(<p)  = one  way  power  antenna  gain 
the  peak; 


normalized  to  one  at 


I 


A0  = clutter  column  width  = 25  mils; 

0.  = azimuth  position  of  a clutter  column  relative  to 
the  antenna  beam  position. 

Notice  that  even  though  g(i)>)  is  nonzero  only  over  16  clutter  columns,  the 
smearing  effect  of  integrating  over  a clutter  column  width  (A0)  for  each 
value  of  0j^,  causes  the  weights,  w(0;j^),  to  have  nonzero  values  over  17  clutter 
columns.  The  number  of  points  in  the  table  equals  the  angular  width  of  17 
columns  divided  by  the  angular  spacing  in  the  table,  i.e.. 


•I 


ki. 


(17) (1/256) 

.-10 


= 68 


The  68  weights  are  written  into  the  random  access  memories  from  the  EROS  digi- 
tal computer  during  simulation  initialization  (immediately  preceding  a simula- 
tion run)  . 

To  avoid  spurious  signals  due  to  azimuth-weight  discontinuities  (refer 
to  Section  2.7),  the  digital  hardware  linearly  interpolates  between  points  in 
the  table*.  This  is,  in  fact,  accomplished  by  interpolating  both  the  azimuth 
readings  and  the  table  values  in  one  operation.  Recall  that  the  azimuth  sam- 
ples are  read  at  a rate  1/256  the  Doppler  sampling  rate  (refer  to  Section  2.9). 
For  each  of  two  successive  azimuth  readings,  the  corresponding  arrays  of  17 
weights. 


Wj,  W2, 


'17 


and 


w. 


w. 


'17 


*Interpolation  is  feasible  for  clutter,  because  only  the  antenna  moves.  It 
would  be  much  more  difficult  for  targets,  where  both  the  antenna  and  the  target 
may  move. 
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are  read  from  the  clutter-azimuth  weight  table.  An  array  weight  increment 


Aw^,  Aw^,  . . . , AWj^^ 


is  then  computed  such  that 


Aw . = 

1 


~ ^i 
256 


The  Interpolated  weights  for  the  17  clutter  columns  can  thus  be  computed 
recursively  for  each  Doppler  cycle  according  to  the  formula 

V 

= w. 


w^(0) 


w^(k  +1)  = w^(k)  + Aw^  j 


i = 1,  . . . , 17 

k = 0,  ....  254 


(3-24) 


3.6  Pseudo  Random  Number  Generator 

The  pseudo  ramdom  number  generator  was  introduced  briefly  in  Section 
3.2,  where  it  was  observed  that  two  points  from  which  to  select  the  random 
inputs  are  sufficient  to  attain  the  desired  amplitude  distribution.  The 
spectral  constraint  on  the  random  number  generator  was  mentioned  in  the 
discussion  of  equation  (3-9);  that  the  autocorrelation  function  have  2 values, 
1 value  at  0-lag  and  a different  value  at  non-O-lag.  There  is  one  additional 
constraint  imposed  by  the  finite  register  length  of  the  filter  implementation: 
that  the  magnitude  Y of  the  random  input  be  sufficiently  small  to  avoid  filter 
overflow.  This  constraint  can  be  expressed  more  precisely*  in  terms  of  the 
filter's  unit-impulse  response  h(nT): 


k < 


Z|h(nT) 


(3-25) 


where  k denotes  the  maximum  magnitude  of  the  total  filter  input.  In  addition 
to  the  pseudo-random  input,  k also  includes  the  noise  due  to  multiplier  trun- 
cation. There  are  2 multipliers  in  the  filter;  if  e^  and  ^2  represent  the 
maximum  magnitudes  assumed  by  their  truncation  errors,  then  (3-25)  becomes 


Y + Of  + 62  < j;[h(nT) 


(3-26) 


The  denominator  in  (3-26)  is  large  for  low  cut-off  frequency  filters,  which 
restricts  y in  this  case  to  small  values  (the  smallest  non-0  value  for  y 
consistent  with  the  16-blt  number  representation  is  2"  ^) . For  higher  cut-off 


*c.f.[6].  Section  9.3.2,  pg.  427. 
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frequency  filters  the  denominator  in  (3-26)  decreases,  and  larger  values  of 
Y can  be  used.  Of  course,  it  would  be  possible  to  use  the  same  "safe"  value 
of  Y (2“^^)  for  high  and  low  frequency  filters.  But  when  the  cut-off  frequency 
is  high,  this  would  reduce  the  filter-output  amplitude  and  restrict  the 
dynamic  range  of  the  simulated  radar  cross-section.  Therefore,  the  pseudo- 
random number  generator  has  been  implemented  with  a controllable  magnitude 
which  can  be  varied  in  powers  of  2.  The  possible  values  of  y range  from 
2“^  through  2“^^. 

The  network  describing  the  pseudo  random  number  generator*  is  illus- 
trated in  Figure  6.  It  consists  of  31  1-bit  delay  elements,  an  exclusive  -OR 
gate,  and  decode  logic  to  produce  the  possible  outputs  +y  or  -y. 


The  generator  has  a cycle  length  of 

,31 

p = 2 - 1 


(3-27) 


I 


[ 


I 


Two  pseudo  random  inputs  are  required  (one  for  each  quadrature  component)  during 
every  Doppler  sampling  cycle  for  each  of  the  544  clutter  cells  in  the  simulated 
antenna  beam  (refer  to  Section  3.5).  The  sampling  period  T]^  of  the  pseudo- 
random number  generator  is  thus  1/1088  of  the  Doppler  sampling  period  T.  A 
portion  of  the  pseudo-random  sequence  used  on  one  of  the  544  filters  would  not 
be  reused  on  any  of  the  other  filters  until  7.9  minutes  later,  and  it  would 
not  be  reused  on  the  same  filter  until  2029  minutes  later.  Therefore,  it  is 
appropriate  to  treat  the  pseudo-random  input  as  aperiodic.  The  relevant  sta- 
tistics of  the  pseudo-random  number  sequence  p(nTj^)  are  listed  as  follows: 

The  mean  of  p(nTj^) 

p{p}  = y(2^^  - 1)  ^ = YP  ^ . (3-28) 

The  variance  of  p(nTj^) 

a^{p}  = Y^  - Y^P~^  . (3-29) 


The  autocorrelation  function  of  p(nTj^) 


Rp(mT^) 


2 n 

Y m = 0 

-Y^p  ^ otherwise  . 


(3-30) 


The  random  Input  to  the  complex  filter  (the  pair  of  identical  real 
filters)  is  denoted  by  r(nT) . The  relevant  statistics  of  r(nT)  are  derived 
from  (3-28),  (3-29),  and  (3-30). 


*The  design  and  the  statistical  properties  of  the  pseudo-random  number  generator 
were  based  on  [7]. 
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EXCLUSIVE-OR 
(ADDITION  MOD-2) 


The  mean  of  r(nT) 


u{r}  = (1  + j)Yp 

The  variance  of  r(nT) 

2r  , 2 2 -2,, 

a {r}  = 2(y  - Y P ) 

The  autocorrelation  function  of  r(nT) 

I 2y^  m = 0 


\(mT) 


-2y^p' 


otherwise 


(3-31) 


(3-32) 


(3-33) 


The  mean  and  variance  of  r enter  into  the  calculations  of  clutter  radar 
cross  section.  The  autocorrelation  function  (3-33)  has  the  requisite  spectral 
property — one  value  at  0-lag  and  a different  value  at  non-O-lag. 

The  evaluation  of  y is  based  on  the  inequality  (3-26).  Expressions  for 
ej^  and  §2  can  be  derived  in  terms  of  the  feedback  multipliers  a,  and  (i2-  The 
sum  |j|h(nT)|  is  also  a function  of  and  02.  However,  a closed  form  expression 
has  not  been  found  for  certain  combinations  of  a-j^  and  02.  Instead  the  sum  is 
evaluated  numerically  using  a software  simulation  of  the  filter.  Section  2.3 
of  the  Eighth  EROS  Quarterly  Report  contains  a more  detailed  discussion  of  this 
calculation. 


3.7  Clutter  Filter  Statistics 


If  the  component  of  clutter  due  to  non-moving  scatterers  is  to  have  a 
specified  radar  cross-section  then  it  is  necessary  to  calculate  the  con- 
stant term  A in  equation  (3-5).  This,  in  turn,  requires  that  the  mean  p of 
the  filter  output  be  known.  Similarly,  equation  (3-6)  indicates  that  in  order 
to  influence  the  radar  cross-section  02  of  the  moving  scatterers  in  clutter, 
it  is  necessary  to  be  able  to  control  the  variance  of  the  filter  outputs. 

Since  this  is  done  by  adjusting  6,  the  post  multiplier  of  the  filter,  an 
expression  relating  6 to  the  variance  of  the  filter  output  is  needed.  The 
development  in  this  section  explains  how  these  statistics  are  obtained. 

The  computations  performed  by  the  clutter  digital  filter  of  Figure  3 
are  described  by  the  difference  in  equations  for  the  states,  Wj  and  W2,  of 
the  two  delay  elements.  The  independent  variable,  time,  is  expressed  in 
multiples  of  the  Doppler  sampling  period  T (refer  to  Section  2.2). 

W^((n  + 1)T)  = a^W^(nT)  + a2W2(nT)  + x(nT)  (3-34) 

W2((n  + 1)T)  = W^(nT)  (3-35) 
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oij^  and  02  are  the  (real)  feedback  multiplier  coefficients.  The  complex  filter 
consists  of  two  identical  real  filters.  Therefore,  it  is  correct  to  treat 
Wj,  W2  , and  X as  complex  numbers.  x(nT)  is  the  sequence  of  random  inputs, 
which  includes  both  the  pseudo-random  numbers  r(nT)  and  the  contributions, 
ei(nT)  and  62 (nT),  due  to  multiplier  truncation  errors. 


x(nT)  = r(nT)  + e^(nT)  + e2(nT).  (3-36) 

The  arithmetic  performed  by  the  digital  filters  employs  16-bit  2's 
complement  binary  fractions;  the  binary  point  is  assumed  to  be  between  the 
sign  bit  and  the  15  fraction  bits.  During  each  Doppler  sampling  cycle,  the 
products  (o]^  - 1)  W-]^  and  02  W2  are  formed,  and  the  30  fraction  bits  of  each 
(real  and  imaginary)  component  of  the  product  is  truncated  to  15  bits;  the  low- 
order  15  bits  are  discarded.  If  [a]  denotes  the  largest  integer  that  does  not 
exceed  a,  and  if  [a  + jb]  denotes  [a]  + jfb],  then  the  truncated  products  are 
expressed  by  t2^^(a]^  - 1)  2“^^  and  [2^^  a2  W2]  2~^^.  Therefore,  the 

truncation  errors  are 

e^  = [2^^  (a^  - 1)  2“^^  - (ct^  - 1)  (3-37) 

and 

®2  ° “2  ^^2^  ~ °'2^2  ’ (3-38) 


The  truncation  error  terms  are  derived  from  the  low-order  bits  of 
product.  In  practice  the  two  multipliers  - 1 and  012  are  different.  This 
implies  that  the  cross  covariance  function  of  e]^(nT)  with  e2(nT)  is  approxi- 
mately 0.  Moreover,  the  multiplicands  Wj^(nT)  and  W2(nT)  generally  change 
values  between  cycles,  which  implies  that  the  autocovariance  functions  of 
ei(nT)  and  e2(nT)  are  approximately  0.  Similarly,  the  cross  covariance  func- 
tions of  e]^(nT)  and  e2(nT)  with  r(nT)  are  approximately  zero.  The  relevant 
statistics  of  x(nT)  follow  from  these  statements  and  from  equations  (3-31), 
(3-32),  (3-33),  and  (3-36). 

The  mean  of  x(nT) 

p{x}  = (1  + j)  YP  ^ + u{ej^}  + p{e2}  • (3-39) 

The  variance  of  x(nT) 

a^{x}  = 2(y^  - Y^  t ^)  + a^{ej^}  + 0^(02)  (3-40) 

The  autocorrelation  function  of  x(nT) 
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m = 0 


(3-41) 


I--  2 2 -2.  ^ 2,  1,2,  1 

2(y  - y p ) + a {e^}  + a {62) 

-2y^  P ^ - 2y^  P ^ + lu{x}l^  otherwise. 

The  expressions  p{e-j^},  p{e2}.  0^(02^],  and  0^(02)  represent  the  mean  and  variance 
of  e^  and  62  respectively.  Equations  (3-39),  (3-40),  and  (3-41)  are  used  in  the 
derivation  of  expressions  for  the  mean  and  variance  of  and  W2.  Moreover, 
equation  (3-41)  expresses  the  requisite  spectral  requirement  for  the  filter 
input — that  its  autocorrelation  function  have  one  value  at  0-lag  and  a different 
value  at  non-O-lag. 

The  mean  and  variance  of  and  W2  are  now  readily  obtained.  Note  that 
W2  is  identical  to  except  for  a delay  of  one  Doppler  sampling  cycle.  There- 
fore, the  statistics  that  describe  W2  are  Identical  to  those  that  describe  Wj^. 
Recall  that  can  be  represented  as  a convolution  of  the  random  process  x 
with  the  unit-impulse  response  h(lT)  of  the  filter.* 

Wj^(nT)  = Z h(iT)  x(nT  - IT)  (3-42) 

The  mean  of  follows  immediately  from  (3-42) 

= F(1  + j)YP~^  + + y{e2}),  (3-43) 


where 


F = Zh(iT). 
i 


(3-44) 


If  G denotes  the  sum  of  squares  of  the  unit-impulse  response 

= Zh^(iT), 


(3-45) 


then  a direct  evaluation  of  the  variance  of  yields 


a^{W^}  = G^a^{x}  + R^(non-0)(F^  - G^)  - F^lwfx}]^ 

= -2F^y^(p“^  + P~^)  + G^(2y^  + a^{e^}  + cr^{e2)  + 2yV^)  • 

The  computations  of  F and  G^  are  simplified  by  the  following  formulas, 

F = .j 

1 - a,  - a- 


(3-46) 


(3-47) 


*Refer  to  Section  8.4  of  [6]. 
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and 


g2  = 


(1  - v 


(1  + 02)1(1  - 02)^  - Oj^^] 


(3-48) 


Equations  (3-43)  and  (3-46)  express  the  mean  and  variance  of  Wj^  as 
functions  of  the  mean  and  variance  of  ej^  and  02-  Conversely,  the  mean  and 
variance  of  ej^  and  e2  can  be  estimated  from  the  mean  and  variance  of  by 
exploiting  the  fact  that  the  quadrature  components  of  Wj^  are  approximately 
Gaussian  distributed.  Therefore,  an  Iterative  algorithm  can  be  employed  to 
determine  the  mean  and  variance  of  Wj^.  The  details  of  this  algorithm  are 
described  in  Section  2.4  of  the  Eighth  EROS  Quarterly  Report.  Included  in 
the  EROS  software  is  a program  to  perform  the  iterations.  Experience  has 
shown  that  convergence  generally  occurs  within  at  most  three  iterations. 


The  weighted  filter  output  Y(nT)  is  related  to  Wj^(nT)  by 

Y(nT)  = 8Wj^(nT) 


(3-49) 


Therefore,  its  mean  and  variance  are  expressed  by 


p{Y}  = 6u{Wj^} 


(3-50) 


and 


a^{Y}  = 


(3-51) 


Combining  (3-51)  with  (3-6)  yields  an  expression  for  8 in  terms  of  02,  the 
radar  cross  section  of  the  moving  clutter  scatterers,  and  the  variance 

of 


0 = 


(3-52) 


(The  constant  K is  defined  in  Section  2.8.) 
o 

The  magnitude  of  the  constant  component  of  clutter  backscatter  is  simply 
/o-j^/Kg  , where  denotes  its  radar  cross  section.  Its  phase  <j)  is  random  (but 
constant),  because  the  constant  component  represents  the  phasor  sum  of  back- 
scatter from  many  reflecting  surfaces  within  the  clutter  cell,  and  each  such 
phasor  is  randomly  oriented.  Thus,  from  equation  (3-5)  an  expression  for  the 
additive  constant  can  be  derived. 


A = e^*  - 0p{W^}  , (3-53) 

where  is  randomly  selected  from  the  uniform  distribution  on  the  interval 

0 < iji  < 2x. 
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4.  TARGET  SIMULATION 


4.1  Target  Data  Management 

The  component  of  radar  backscatter  produced  by  targets  is,  as  for  clutter, 
represented  by  a complex  envelope  in  digital  form.  The  Doppler  sampling  fre- 
quency, 4.167  KHz,  permits  the  unambiguous  representation  of  radial  velocities 
up  to  approximately  30  meters/second  or  67  miles/hour.* 

The  samples  of  the  target  signal  are  obtained  from  analog  recordings  of 
coherently  detected  backscatter.  Each  sample  consists  of  an  in-phase  and  a 
quadrature  component.  For  this  purpose  EROS  contains  two  A/D  converters,  asso- 
ciated sampling  and  synchronization  hardware,  and  software  to  store  the  samples 
on  computer  disk.  For  scaling  purposes  the  user  specifies  the  radar  cross 
section  of  the  recorded  target. 

In  addition  to  sampled  backscatter,  target  simulation  also  requires  a 
description  of  maneuvers — the  path  on  which  the  target  moves  and  an  association 
of  points  on  the  path  with  times  during  the  simulation.  The  path  is  segmented 
into  "legs";  the  two  endpoints  of  each  leg  are  called  "nodes".  The  target 
motion  while  traversing  a leg  is  constant  radial  velocity  and  constant  angular 
velocity.  Legs  form  straight  line  segments  on  a B-scope  display.  For  each 
leg  the  following  information  is  provided: 

° the  (range,  azimuth)  coordinates  of  the  nodes; 

° the  simulation  times  associated  with  the  nodes; 

o the  name  of  a target  recording  or  the  specification  that 
the  target  has  disappeared  during  this  leg; 

o the  portion  of  the  recording  to  use; 

o a weight  (amplification  or  attenuation  factor)  to  be 
applied  to  the  radar  cross  section  for  the  duration 
of  the  leg. 

The  weight  can  be  used  to  account  for  partial  masking. 

A computer  program  called  "target  compilation"  combines  this  data  and 
stores  it  on  disk  in  a form  that  is  optimal  for  real-time  retrieval.  For  each 
Doppler  sampling  cycle  the  target's  azimuth  and  range  are  computed.  Seventeen- 
bit  azimuth  accuracy  is  retained  in  accordance  with  Section  2.7.  The  instant- 
aneous range  of  the  target  is  expressed  in  terms  of  range-ring  numbers;  the 
user  will  have  specified  the  range  from  the  radar  to  the  closest  simulated 


*This  is  based  on  the  AN/PPS-15  radar's  transmitter  frequency  of  10.2  to 
10.4  GHz. 


ring,  and  the  rings  are  15.027  meters  wide  (see  Section  2.1).  The  apportion- 
ment of  the  target  sample  between  two  adjacent  range  rings  has  been  described 
in  Section  1.2.  Each  Doppler  sampling  cycle  is  represented  by  four  12-blt 
sample  components:  the  portions  of  the  in-phase  and  quadrature  components 
in  range-ring  r and  in  range-ring  r + 1.  These  samples  are  stored  in  four 
16-bit  words.  The  remaining  16  bits  not  occupied  by  the  sample  components 
contain  range-ring  data,  azimuth  data,  and  synchronization  information. 

During  real-time  simulation  the  computer  retrieves  the  target  data, 
compares  its  azimuth  with  that  of  the  antenna  beam,  obtains  an  azimuth  weight 
from  an  antenna-pattern  table,  and  transmits  information  to  the  digital  hard- 


® the  antenna  pattern  weight; 
o the  in-phase  and  quadrature  components; 

° the  range-ring  number; 

° synchronization  information. 

The  digital  hardware  multiplies  the  sample  components  by  the  weight  and 
combines  the  products  with  the  clutter  samples  at  the  appropriate  range  rings. 


4.2  Real-Time  Limitations 

The  time  interval  of  240  usee  between  target  samples  imposes  a signifi- 
cant burden  on  the  computer  during  real-time  simulation.  Much  of  the  effort 
in  preparing  a target  sample  is  expended  in  interacting  with  the  subject 
radar's  antenna.  It  would  be  impossible  for  the  computer  to  perform  the 
antenna-pattern  multiplications  in  real-time.  For  this  reason  the  target 
azimuth  weighting  is  done  in  the  special  purpose  digital  hardware. 

Even  without  having  to  perform  the  multiplications,  a substantial 
amount  of  time  is  required  in  each  Doppler  sampling  cycle  to  look  up  the 
azimuth  weight  and  to  transmit  the  data  to  the  digital  hardware.  Timing 
measurements  indicate  that  more  than  100  psec  are  needed  to  perform  the  com- 
putation for  one  target  during  each  240  ysec  Doppler  cycle.  Since  the  computer 
has  to  provide  clutter  parameters  in  addition  to  the  target  samples,  the  com- 
puter would  not  be  capable  of  simulating  more  than  one  target  in  the  antenna 
beam  at  a time. 

The  following  restrictions  have  been  imposed  on  the  target  portion  of 
the  scenario  to  meet  the  real-time  requirements: 

o A target  scenario  can  consist  of  one  or  two  target  paths. 

The  target  paths  consist  in  general  of  several  legs. 

Associated  with  each  target  path  are  sectors  defined  by  a 
beginning  azimuth  and  an  ending  azimuth.  All  nodes  of 
the  path  must  fall  between  these  two  azimuths. 


o Each  target  path  can  contain  several  targets,  but  only  one 
target  at  a time.  The  legs  are  simulated  in  sequence. 

" The  angle  between  the  beginning  and  ending  azimuth  must  be 
2734  mils  or  less. 

° The  sectors  of  the  two  paths  must  not  overlap.  The  sectors 
must  be  separated  by  at  least  500  mils. 


4.3  Target  Sample  Blocking 

Consider  the  schedule  for  target  processing  shown  in  Figure  7(a).  Real- 
time processing  is  subdivided  into  240  ysec  intervals.  The  target  data  for 
each  Doppler  sampling  cycle  are  sent  to  the  digital  hardware  at  the  beginning 
of  each  of  these  intervals.  The  remaining  time  at  the  end  of  the  240  Msec 
interval  is  available  for  supplying  clutter  parameters.  In  order  to  alert 
the  computer  to  start  target  processing  at  the  beginning  of  each  interval, 
either  Interrupts  must  be  provided  at  a 4.167  KHz  frequency,  or  the  clutter 
program  must  check  a status  flag  frequently.  The  interrupt  method  is  unsatis- 
factory, because  the  time  required  to  switch  status  (including  saving  and 
restoring  contents  of  six  general-purpose  registers)  would  severely  aggravate 
the  real-time  load.  The  frequent  checking  of  the  status  flag  within  the  short 
time  available  at  the  end  of  the  240  ysec  interval  is  similarly  infeasible. 

An  alternative  schedule  which  resolves  these  difficulties  is  shown  in 
Figure  7(b).  Real-time  processing  is  subdivided  into  Intervals  of  length 
240  n Msec,  where  n equals  128  for  the  current  version  of  EROS.  The  target 
data  for  n Doppler  sampling  cycles  are  sent  to  the  digital  hardware  at  the 
beginning  of  each  of  these  intervals.  The  remaining  time  at  the  end  of  each 
interval  is  now  in  a more  usable  form  for  performing  clutter  processing.  The 
digital  hardware  alerts  the  computer  at  the  end  of  each  240  n Msec  interval 
by  setting  a status  bit. 

Buffer  storage  must  be  provided  to  hold  a block  of  data  since  it  is 
generated  well  in  advance  of  the  time  the  radar  would  receive  that  target  sig- 
nal. The  storage  could  have  been  part  of  the  computer’s  core  memory  instead 
of  dedicated  random-access  memory  in  the  digital  hardware.  The  latter  choice 
was  made  for  a number  of  reasons,  including  software  simplicity,  elimination 
of  need  for  direct-memory-access  hardware,  and  reduction  of  the  computer's  core- 
memory usage.  Two  such  buffers  are  provided,  so  that  the  computer  can  be 
writing  into  one  of  them  while  the  hardware  reads  previously  stored  data  from 
the  other  one. 
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TARGET  NON-TARGET 


5. 


EROS  IMPLEMENTATION 


This  section  presents  brief  descriptions  of  the  major  units  comprising 
the  EROS  implementation. 

5.1  Simulation  Preparation  Software 

Simulation  preparation  software  consists  of  programs  to  assist  the  oper- 
ator in  describing  the  scenario.  Much  of  the  information  is  maintained  in 
reference  files  (called  "libraries”) , whose  contents  may  have  widespread  appli- 
cability in  many  scenarios.  All  input  entered  by  the  user  is  checked  for  valid- 
ity, and  diagnostics  are  written  when  errors  are  detected.  Generous  provision 
is  made  for  the  user  to  append  comments  to  document  his  simulation  data. 

Optional  printouts  permit  the  user  to  generate  reference  reports  that  define 
the  files'  contents  and  that  contain  the  user-written  comments. 

The  records  in  the  clutter  spectral  library  specify  the  spectral  prop- 
erties C'f  the  clutter  filters.  Each  record  corresponds  to  a different  spectral 
behavior,  and  is  identified  by  a unique  integer.  Facilities  are  provided  for 
the  EROS  user  to  add  or  change  records  in  this  library. 

The  complete  clutter  library  consists  of  records  which  contain  both 
spectral  and  radar  cross  section  information.  To  define  a record  in  this 
library,  the  user  supplies  (1)  a 2-character  code  to  identify  the  record; 

(2)  an  integer  that  refers  to  a record  in  the  clutter  spectral  library;  and, 

(3)  radar  cross  section  values  for  the  fixed  and  moving  components  of  clutter. 
The  clutter  library  update  program  retrieves  the  spectral  information  inferred 
by  the  given  spectral-reference  Integer.  Facilities  are  provided  for  the  EROS 
user  to  add,  delete,  or  change  records  in  the  library. 

The  clutter  scenario  description  program  assists  the  user  in  defining 
the  positional  dependence  of  clutter.  The  computer's  graphic  cathode-ray-tube 
and  light  pen  are  made  available  to  the  user  to  expedite  his  definition  of  the 
8192-cell  clutter  description  array.  The  type  of  clutter  in  each  cell  is 
identified  by  the  two-character  code  that  refers  to  a record  in  the  complete 
clutter  library.  The  user  may  also  need  to  define  a clutter  weight  array, 
which  specifies  an  amplification  or  attenuation  to  be  applied  to  the  radar  cross 
section  in  each  cell.  The  weights  can  account  for  partial  masking,  range 
dependence  of  the  cell's  areas,  incidence  angle  effects,  etc.  A unique  file 
name  is  associated  with  each  clutter  description  array  and  each  clutter  weight 
array.  Facilities  are  provided  for  the  EROS  user  to  revise  an  array  or  to 
define  a new  array  as  a modification  to  a previously  defined  array.  The  clutter 
scenario  compilation  program  combines  specified  arrays  and  reference-file 
information  to  produce  the  "compiled  clutter  scenario  file"  in  an  optimal  format 
for  real-time  simulation. 


The  target  portion  of  a scenario  is  built  up  from  recordings  and  from 
target  path  definitions.  The  recorded  information  extracted  from  the  A/D 
converters  is  stored  on  disk  files  as  described  in  Section  4.1.  A unique  file 
name  is  associated  with  each  recording.  The  computer's  cathode-ray-tube  and 
light  pen  are  available  to  the  user  for  describing  the  target  paths,  which 
are  depicted  in  the  form  of  a B-scope  display.  The  target  compilation  program 
combines  one  or  two  target  paths  and  the  recordings  referenced  in  these  paths. 
It  produces  the  "compiled  target  scenario  file"  in  an  optimal  format  for 
real-time  simulation. 

Antenna  patterns  are  represented  during  real-time  simulation  by  a 
4096-point  table  for  attenuating  target  signals  and  by  a 68-point  table  for 
attenuating  clutter  signals.  The  user  specifies  the  pattern  with  13  points 
spaced  at  1-degree  intervals  from  the  beam  center.  The  computer  translates 
this  specification  into  an  antenna  pattern  file  which  holds  the  two  tables; 
formula  (3-23)  is  used  to  calculate  the  clutter  antenna  pattern  table. 

The  range-weight  file  consists  of  197  attenuation  factors  corresponding 
to  range  rings  from  60  to  3000  meters  at  15-meter  intervals.  Specified  range- 
weight  Inputs  are  converted  to  fixed-point  format  for  real-time  simulation. 
Note  that  only  64  of  these  weights  are  used  during  any  given  simulation  run. 
The  minimum  range  for  the  simulation  determines  which  of  these  64  weights  to 
use. 


5.2  Real-Time  Software 


The  user  initiates  a simulation  run  by  entering  the  name  of  an  "PROS 
scenario  file".  EROS  scenario  files,  which  are  created  by  the  user  via  the 
operating  system's  text  editor,  contain  the  names  of  the  component  files: 

^ a compiled  clutter  scenario  file, 

o a compiled  target  scenario  file, 

o an  antenna  pattern  file, 

o a range-weight  file. 

During  simulation  initialization,  various  tables  and  constants  are  read  into 
the  computer  and  into  digital-hardware  memories.  These  include  the  antenna- 
pattern  and  range-weight  tables.  EROS  starts  producing  simulated  backscatter 
when  the  radar's  antenna  reaches  a specified  azimuth  and  is  proceeding  in  a 
specified  scan  direction.  This  increases  the  repeatability  of  EROS  experiments. 
After  simulation  has  started,  the  real-time  program  responds  to  the  antenna's 
azimuth  by  copying  appropriate  clutter  and  target  data  from  the  disk  files 
to  the  digital  hardware. 


5.3  Digital  Hardware 


The  digital  hardware  performs  those  functions,  where  the  speed  and  the 
number  of  computations  are  beyond  the  real-time  capability  of  the  computer. 
The  data  Inputs  to  the  digital  hardware  are  supplied  by  the  computer: 

o filter  parameters  and  current  azimuth  for  clutter; 

o signal  samples,  range,  and  azimuth  weight  for  targets. 

The  clutter  processor  continuously  calculates  synthetic  clutter  back- 
scatter  from  the  544  cells  in  the  currently  simulated  antenna  beam.  Each 
Doppler  sampling  cycle  consists  of  performing  one  computation  step  for  each 
filter,  applying  the  antenna  pattern  weights,  and  performing  azimuth  inte- 
gration. The  azimuth  input  from  the  computer  is  used  to  determine  which 
antenna-pattern  weight  applies  to  a given  cell.  Azimuth  integration  consists 
of  adding  together  the  returns  from  cells  at  the  same  range.  The  result  is 
a sampled  profile  of  clutter  return 

Sq.  . . . s^3 

spaced  at  15-meter  intervals,  such  that 
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The  values  of  these  64  samples  are  recalculated  each  240  psec. 

The  target  processor  applies  the  azimuth  weight  to  the  two  quadrature 
components  of  each  target  sample  and  stores  the  products  with  the  associated 
range-ring  number  in  a random  access  memory  (RAM).  In  general,  two  such  sets 
of  data  are  required  for  each  Doppler  sampling  cycle,  because  part  of  a target 
sample  is  allocated  to  one  range  ring,  and  the  balance  is  allocated  to  an 
adjacent  range  ring.  The  RAM  is  capable  of  storing  data  for  128  Doppler  samp- 
ling cycles.  Status  bits  are  stored  in  the  RAM  to  denote  the  end  of  each 
Doppler  sampling  cycle  and  the  end  of  128  Doppler  sampling  cycles.  Two  such 
RAM's  are  used,  so  that  the  computer  can  load  one  of  them  while  the  target 
processor  retrieves  previously  loaded  data  from  the  other  one.  Each  retrieved 
sample  is  added  to  the  clutter  range  profile  at  the  range  specified  by  the 
retrieved  range-ring  number. 

The  video  signal  is  produced  by  storing  each  range  profile  in  a RAM 
and  by  reading  it  back  at  the  radar's  pulse  repetition  frequency  (PRF) . Since 
the  Doppler  sampling  frequency  is  1/9  the  PRF  of  the  AN/PPS-15  radar  (see 
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Section  2.2),  each  range  profile  is  read  back  nine  times.  Two  RAM's  are 
used  for  this  purpose,  so  that  one  of  them  can  be  loaded  from  the  target  pro- 
cessor while  the  video  signal  is  read  from  the  other  one. 

Four  clocks  are  used  to  synchronize  the  digital  hardware.  The  clutter 
and  target  processors  employ  a 10  MHz  clock  derived  from  a crystal  oscillator 
and  a 4.167  KHz  clock  derived  by  dividing  the  radar's  37.5  KHz  PRF  by  9.  The 
video  signal  generator  employs  the  4.167  KHz  clock  for  loading  the  RAM's;  for 
reading  back  the  video  signal  it  employs  the  37.5  KHz  PRF  clock  and  a 10  MHz 
clock  derived  from  a quenched  oscillator.  The  10  MHz  frequency  corresponds 
to  the  15-meter  range-ring  spacing.  The  quenched  oscillator  is  restarted  by 
a signal  derived  from  the  radar's  PRF  trigger  and  provides  a stable  repre- 
sentation of  the  delay  between  radar  transmission  and  reception. 


5.4  Analog  Hardware 

Analog  hardware  is  needed  in  any  radar  signal  simulation  to  convert 
digital  samples  into  a form  suitable  for  insertion  into  the  radar.  Moreover, 
in  EROS  the  range-weighting  and  range-integration  functions  are  performed 
in  analog  hardware  to  meet  the  AN/PPS-15  radar's  dynamic-range  and  video- 
bandwidth capabilities.  A preliminary  design*,  in  which  these  computations 
were  to  be  performed  digitally,  was  rejected  in  order  to  improve  signal 
quality  and  to  simplify  hardware  and  software. 

Twelve-bit  D/A  converters  are  available  with  speeds  adequate  to  process 
real-time  video  signals.  Three  such  converters  are  used;  (1)  to  convert 
sign  and  magnitude  of  In-phase  (I)  data  prior  to  range  weighting;  (2)  to 
convert  sign  and  magnitude  of  Quadrature  (Q)  data  prior  to  range  weighting; 
and  (3)  to  convert  magnitude  of  the  range  weight  values.  The  analog  output 
of  the  range  weight  D/A  converter  is  multiplied  by  the  I signal  and  by  the 
Q signal  in  separate  analog  multipliers. 

For  video  signal  simulation,  digital  signals  representing  the  I and  Q 
content  and  range  weight  of  successive  range  rings  are  latched  into  the  three 
D/A  converters  at  100  nanosecond  intervals.  Outputs  of  the  analog  multipliers 
are  integrated  in  analog  integrators,  which  are  reset  periodically  by  pulses 
coincident  with  the  transmitter  modulation.  Since  signals  from  all  range 
rings  are  present  for  equal  time  increments,  integration  with  respect  to  time 
is  equivalent  to  algebraic  addition  of  the  data  from  successive  range  rings. 
The  output  from  the  integrators  is  a faithful  representation  of  the  envelope 
of  AN/PPS-15  radar  return,  because  the  radar  transmits  CW  with  180°  phase 
reversals  at  the  PRF.  The  integrators  are  also  quite  effective  in  reducing 
the  bit-switching  transients  that  are  inherent  in  D/A  converters.  The  output 
of  the  Integrators  are  fed  into  variable-gain  video  amplifiers,  which  provide 
calibration  and  the  addition  of  adjustable  simulated  receiver  noise.  Simple 
band-pass  filters  are  used  to  limit  the  video  signal  bandwidth  to  that  of  the 
radar  pre-amplifier,  and  to  further  attenuate  bit-switching  transients  that 
originate  in  the  D/A  converters. 


*Refer  to  Section  4 of  the  Fourth  EROS  Quarterly  Report. 
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For  simulation  of  range-gated  Doppler  signals,  the  EROS  operator  selects 
one  of  the  64  range  rings  by  means  of  a knob  on  the  analog  hardware.  Digital 
signals  representing  the  1 content,  Q content,  and  range  weight  of  the  selected 
range  ring  are  latched  into  the  D/A  converters  at  the  transmitter  modulation 
rate,  and  each  sample  is  held  for  a time  equal  to  the  interval  between  trans- 
mitter modulations.  The  range  Integrators  have  been  switched  out,  and  the 
range  weighting  analog  multipliers  are  connected  through  low-pass  filters  to 
variable-gain  Doppler  output  amplifiers. 

The  present  method  of  analog  range  weighting  lends  itself  readily  to 
simulation  of  other  types  of  radar  modulation.  In  particular,  the  bipolar 
video  signal  of  conventional  pulse-Doppler  radar  differs  from  the  AN/PPS-15 
video  in  that  the  returns  from  a single  range  resolution  cell  contribute  to 
the  video  signal  at  any  given  point  in  the  Interpulse  interval.  This  signal 
is  simulated  in  EROS  simply  by  switching  out  the  range  integrator  and  adjusting 
signal  and  noise  levels  appropriately  for  the  parameters  of  the  radar  under 
test.  The  range-gated  Doppler  signal  for  pulse  modulation  and  for  the 
AN/PPS-15  modulation  do  not  differ  in  form,  and  only  the  adjustment  of  signal 
and  noise  levels  is  required. 


5.5  Current  Status 


Each  of  the  component  units  of  EROS  has  been  fabricated  and  unit 
tested.  The  software  has  been  written,  debugged,  and  system  tested.  However, 
the  integration  of  these  components  is  still  in  progress  and  a complete  system 
test  will  be  conducted  under  the  Contract  DAAB07-77-C-2174 , which  is  an 
extension  of  the  current  program. 
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EFFECT  OF  DISCRETE  REPRESENTATION  OF  ANTENNA  PATTERN 
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The  effect  of  the  antenna  pattern  in  the  horizontal  plane  is  simulated 
in  EROS  by  applying  an  appropriate  weight  to  the  radar  return  from  each  back- 
scatter  element.  The  weight  depends  upon  the  difference  between  the  azimuth 
angle  describing  the  location  of  the  backscatter  element  and  the  azimuth  angle 
of  the  antenna  direction.  The  value  of  the  weight  equals  the  two-way  voltage 
antenna  gain  function  evaluated  at  this  angle  difference.  A tabular  repre- 
sentation of  the  antenna  gain  function  will  be  used  in  EROS  for  two  reasons: 

(1)  the  computer  can  perform  table  lookups  faster  than 
calculate  a formula; 

(2)  greater  flexibility  is  permitted  in  representing 
the  antenna  pattern;  e.g,,  measured  antenna  gain 
data  can  be  applied  to  the  simulation  even  if  the 
data  cannot  be  conveniently  reduced  to  a formula. 

A small  number  of  table  entries  (e.g.,  less  than  50)  would  be  adequate 
for  accurate  representation  of  antenna  gain  if  linear  interpolation  were  used. 
However,  the  computation  time  required  to  perform  even  linear  interpolation 
would  add  a significant  burden  to  the  computer’s  real-time  load.  To  avoid 
interpolation,  we  would  design  the  antenna  pattern  table  with  a large  number 
of  values  of  the  gain  function  at  evenly  spaced  angles.  The  real-time  target- 
processing program  would  then  select  the  table  value  whose  corresponding  angle 
is  closest  to  the  target's  azimuth  difference. 

Consider  the  effect  of  scanning  at  constant  angular  rate  across  a 
stationary  reflector  with  large  radar  cross  section  (e.g.,  a corner  reflector). 
The  radar  return  signal  is  amplitude  modulated  by  the  antenna  gain  function, 
while  the  angle  is  varied  at  the  antenna  scan  rate.  Let  a denote  the  angular 
rate  of  antenna  scan,  and  A9  denote  the  angular  separation  between  adjacent 
points  in  the  antenna  gain  function  table.  Using  this  tabular  representation 
of  antenna  gain,  an  EROS  simulation  of  the  corner  reflector  scenario  would 
produce  the  staircase  function  shown  in  Figure  8.  The  simulated  signal  remains 
constant  for  a duration  of  and  then  jumps  to  the  next  value  in  the  table. 

Of  course,  if  ^ is  less  than  the  240  usee  sampling  interval,  the  simulation 
is  perfectly  adequate.  However,  for  larger  values  of  &9  the  jump  discontin- 
uities increase  in  magnitude  and  decrease  in  frequency,  and  they  introduce 
spurious  signals  at  frequencies  within  the  pass  band  of  the  AN/PPS-15.  It  is 
therefore  of  Interest  to  determine  the  largest  tolerable  value  of  A6  for  which 
the  spurious  signal  is  at  a svif f iciently  low  power.  The  following  analysis 
addresses  this  question. 


In  performing  this  analysis  several  assumptions  were  made: 

® The  corner  reflector  scenario  is  the  worst-case  problem. 

In  other  words,  the  spurious  signal  introduced  by  the 
discrete  antenna  pattern  representation  with  other  illu- 
minated targets  will  be  tolerable  if  the  corner  reflector 
problem  is  solved. 

» If  the  power  of  the  spurious  signal  is  at  least  72  dB 
lower  than  the  saturating  signal  level,  it  has  been 
adequately  suppressed. 

o Frequency  components  of  the  spurious  signal  above  3 kHz 
can  be  ignored,  because  they  will  be  filtered  out  by 
the  digital  interpolation  up  to  37.5  kHz,  the  post-audio 
filter  in  the  EROS  analog  hardware,  and  audio  filtering 
within  the  AN/PPS-15  radar. 

» The  antenna  scans  at  a fixed  angular  rate  a equal  to 
90  mils/second. 

o The  largest  value  of  the  slope  of  the  normalized  antenna 
pattern  (whose  gain  varies  from  0 at  the  nulls  to  1 at 
the  nose)  is  .0138  This  value  was  derived 

from  a sample  AN/PPS-15  measured  antenna  pattern. 

The  spurious  signal  introduced  by  the  discrete  antenna  gain  function  is 
the  difference  between  the  staircase  function  and  scan  modulation  function  in 
Figure  8.  This  spurious  signal  has  the  shape  of  a sawtooth  wave  with  amplitude 
varying  as  the  slope  of  the  scan  modulation  function.  Clearly  the  power  in 
the  sawtooth  signal  is  maximum  at  the  steepest  slope.  Therefore,  we  simplify 
the  analysis  by  considering  the  periodic  sawtooth  function  E(t)  illustrated 
in  Figure  9.  whose  amplitude  a equals  1/2  the  maximum  discontinuity  in  the 
staircase  function. 

Let  W(0)  denote  the  two-way  voltage  antenna  gain  function  of  azimuth 
angle  9,  and  let  A denote  the  amplitude  of  the  radar  return  when  the  antenna 
is  pointing  directly  at  the  corner  reflector.  The  magnitude  of  the  jump  dis- 
continuities in  the  staircase  function  of  Figure  8.  are  thus  equal  to 
a|w(9  + A0)  - W(0)|.  Since  EROS  scenarios  will  not  include  reflectors  whose 
radar  return  exceeds  the  saturating  signal  level,  we  can  safely  assume  that 
A does  not  exceed  the  saturating  signal  level  • Therefore, 

a - 1 A|W<9  40)  - f A X 40  < j X 40  . (Al) 
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The  Fourier  series  representation  of  E(t)  is 


where 


E(t) 


na-co 


f 

o 


g 

A0 


(A2) 


Note  that  f denotes  the  frequency  at  which  the  antenna  gain  function  is  sampled. 
After  filtering  out  components  above  3 kHz,  the  low-pass  filtered  error  func- 
tion [E(t)]Lp  is  expressed  by 


[E(t)  ]j^p 


(-1) 


n ai  gj  a 2Tr 

niT 


where  k is  largest  integer  satisfying 

f k < 3000  . (A3) 
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By  Parseval's  theorem  the  power  P in  [E(t)]p^p  is  given  by 
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n^O 
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Since  the  power  in  the  saturating  signal  is  A^^^  • since  P must  be  at 

least  72  dB  below  the  power  of  the  saturating  signal, 
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< -72 


Combining  (Ai),  (A2) , (A4) , and  (A5) , we  get 
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< -72 


(A5) 


(A6) 


Finally  by  substituting  the  values  of  X and  a stated  in  the  assumptions,  we 
compute  the  smallest  f^  satisfying  (A3)  and  (A6)  to  be  1244  Hz.  By  (A2)  the 
corresponding  angular-resolution  constraint  is  that  A9  < .0723  mils  = 1.13  x 
2^q“5  circles.  Since  angles  will  be  expressed  during  EROS  simulation  in  terms 
of  binary  fractions  of  a circle,  at  least  17  bits  are  required  to  attain  this 
angular  resolution. 
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